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By confining water in nano-pores of silica glass, we can bypass the crystallization and study the 
pressure effect on the dynamical behavior in deeply supercooled state using neutron scattering. We 
observe a clear evidence of a cusp-like fragile-to-strong (F-S) dynamic transition. Here we show 
that the transition temperature decreases steadily with an increasing pressure, until it intersects the 
homogenous nucleation temperature line of bulk water at a pressure of 1600 bar. Above this pressure, 
it is no longer possible to discern the characteristic feature of the F-S transition. Identification of 
this end point with the possible second critical point is discussed. 
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Water is a continuing source of fascination to scientists 
because of its many counterintuitive low-temperature 
properties. Although the stable form of water at suffi- 
ciently low temperature is inevitably crystalline, liquid 
water can also exist in a metastable supercooled state 
far below the melting point. One of the most intriguing 
questions related to the unusual properties of supercooled 
water is whether two critical points may exist in a single 
component liquid 0. 

It has been known since 1970s that a number of 
thermodynamic response functions of supercooled wa- 
ter, notably the isothermal compressibility Kt and the 
constant-pressure specific heat Cp, show a power-law 
divergence behavior at a singular temperature, experi- 
mentally determined to be 228 K at the ambient pres- 
sure Concurrently, the transport properties, such 
as the shear viscosity r\ and the inverse self-diffusion con- 
stant D, diverge according to power-laws toward the same 
singular temperature 0,0 • The anomalies of the thermo- 
dynamic quantities become plausible if one postulates the 
existence of a second low-temperature critical point at 
about 228 K and at somewhat elevated pressure ■ On 
the other hand, the transport coefficient anomalies are 
reminiscent of the dynamical behavior of a supercooled 
liquid near the so-called kinetic glass transition temper- 
ature, predicted by Mode-Coupling theory (MCT) jg. 

Search for the predicted first-order liquid-liquid 
transition line and its end point, the second low- 
temperature critical point 00 in water, has been ham- 
pered by intervention of the homogenous nucleation pro- 
cess, which takes place at 235 K at the ambient pressure. 
However, by confining water in nano-pores of mesoporous 
silica MCM-41-S with cylindrical pores of 14 A diame- 
ter, we have been able to study the dynamical behavior 
of water in a temperature range down to 160 K, with- 
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out crystallization. Using high-resolution Quasi-Elastic 
Neutron Scattering (QENS) method and Relaxing-Cage 
Model (RCM) for the analysis, we determine the tem- 
perature and pressure dependences of the average trans- 
lational relaxation time (tt) for the confined supercooled 
water. 

Micellar templated mesoporous silica matrices MCM- 
41-S, which have 1-D cylindrical pores arranged in 2-D 
hexagonal arrays, were synthesized by following a sim- 
ilar method for synthesizing MCM-48-S previously 0. 
Using the short chain cationic surfactants, C12TMAB, 
as tem plat es and adding /3-type zeolite seeds as silica 
source [l(J, we obtain MCM-41-S (S denotes seed) with 
smaller pore sizes and stronger silica walls than the tra- 
ditional MCM-41. In this pressure experiment, we chose 
the silica with a pore diameter of 14 A because for the 
fully hydrated sample the differential scanning calorime- 
try (DSC) data shows no freezing peak down to 160 K. 
The sample is then hydrated by exposing to water vapor 
in a closed chamber until it reaches the full hydration 
level of 0.5 gram H2O/I gram Silica. 

High-resolution QENS spectroscopy method is used to 
determine the temperature and pressure dependences of 
(tt) for the confined water. Because neutrons can easily 
penetrate the thick-wall high-pressure cell and because 
it is predominantly scattered by hydrogen atoms in wa- 
ter, rather than by the matrices containing it, incoher- 
ent QENS is an appropriate tool for our study. Using 
two separate high-resolution QENS spectrometers, we are 
able to measure the translational -relaxation time from 
0.2 ps to 10,000 ps over the temperature and pressure 
range. The high-pressure experiments were performed 
at both the High-Flux Backscattering (HFBS) and the 
Disc-Chopper Time-of-Flight (DCS) spectrometers in the 
NIST Center for Neutron Research (NIST NCNR). The 
two spectrometers used to measure the spectra have two 
widely different dynamic ranges (for the chosen experi- 
mental setup), one with an energy resolution of 0.8 fieV 
(HFBS) and a dynamic range of ±11 iieV [ll|, 

and the 

other with an energy resolution of 20 fieV (DCS) and a 
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dynamic range of ±0.5 meV [l2( in order to be able to 
extract the broad range of relaxation times from the mea- 
sured spectra. The same high pressure system, including 
specially designed aluminium pressure cell, was used on 
the two instruments. Helium gas, the pressure-supplying 
medium, fills the whole sample cell, and applies pressure 
to the fully hydrated sample. The experiment at each 
pressure was done with a series of temperatures, cover- 
ing both below and above the transition temperature. 
Altogether, 1000 spectra were collected, spanning 9 pres- 
sures: ambient, 100, 200, 400, 800, 1200, 1600, 2000, and 
2400 bars. 

QENS experiments measure the Fourier transform of 
the Intermediate Scattering Function (ISF) of the hy- 
drogen atoms, FH{Q,t), of water molecules in the pores 
of the silica matrix. Molecular Dynamics (MD) simu- 
lations have shown that the ISF of both bulk [13j and 
confined [14| supercooled water can be accurately de- 
scribed as a two-step relaxation: a short-time Gaussian- 
like (in-cage vibrational) relaxation followed by a plateau 
and then a long-time (time > 1.0 ps) stretched expo- 
nential relaxation of the cage. The so-called Relaxing 
Cage Model (RCM) 0], which we use for data analy- 
sis, models closely this two-step relaxation and has been 
tested extensively against bulk and confined supercooled 
water through MD and experimental data 0, 0, ^( . 
By considering only the spectra with wave vector trans- 
fer Q < 1.1 A" 1 , we can safely neglect the contribution 
from the rotational motion of water molecule in ISF |l5j | . 
The RCM describes the translational dynamics of water 
at supercooled temperature in terms of the product of 
two functions: 



F H (Q,t) « F T (Q,t) = F s (Q,t)exp[-(t/T T (Q)) 
tt(Q) = T (0.5QH , (r T ) =r r(l//3)//3, 



(1) 



where the first factor, F s (Q,t), represents the short- 
time vibrational dynamics of the water molecule in the 
cage. This function is fairly insensitive to temperature 
variation, and thus can be calculated from MD simula- 
tion. The second factor, the a-relaxation term, contains 
the stretch exponent /3, and the Q-dependent transla- 
tional relaxation time tt (Q), which is a strong function 
of temperature. The latter quantity is further specified 
by two phcnomcnological parameters tq and 7, the expo- 
nent controlling the power-law Q-dependence of tt (Q)- 
(tt) is a Q-independent quantity where T is the gamma 
function. It essentially gives a measure of the structural 
relaxation time of the hydrogen-bond cage surrounding a 
typical water molecule. The temperature dependence of 
the translational relaxation time at each pressure is then 
calculated from three fitted parameters, To, /3, and 7, by 
analyzing a group of quasi-elastic peaks at different Q 
values simultaneously. For this analysis, we chose seven 
spectra for data taken at HFBS and eleven spectra for 
data taken at DCS, at each temperature. 

We show in Fig^ as an example, two complete sets 
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FIG. 1: Figs[TJ\ and [Tpl (left panels) show QENS spectra 
measured at Q — 0.58 A -1 , at two pressures, 800 bar and 
1600 bar, and at a series of temperatures. Figs fit? and 
(right panels) show the RCM analysis of one of the spectrum 
from each pressure. The resolution function in each case is 
shown by a dashed line. 



(temperature series) of QENS area-normalized spectra. 
The broadening of the quasi-elastic peaks becomes more 
and more noticeable as temperature increases. In Fig^\, 
we may notice, from shoulders of these spectral lines, that 
two groups of curves, 231-250 K and 198-209 K, are sepa- 
rated by the curve at a temperature of 220 K. This visual 
information reinforces the result of the detail analysis 
shown in Fig|2 that there is an abrupt dynamical tran- 
sition at T L = 216 K. Fig[Tp shows the RCM analysis of 
the spectrum taken at T — 220 K, close to the transition 
temperature. On the other hand in FigQp, the spectra 
at pressure 1600 bar show a rather smooth variation with 
temperature, indicating that there is no sharp transition. 
This pressure corresponds to the end point of the line of 
F-S transitions shown in Fig|21 Fig^) is, again, an RCM 
analysis of the spectrum taken at T = 214 K of this pres- 
sure. RCM, as one can see, reproduces the experimental 
spectral line shapes of confined water well. The broaden- 
ing of the experimental data over the resolution function, 
shown in Figs ^3 and Fig^D, leaves enough dynamic in- 
formation to be extracted by RCM. 

The behavior of shear viscosity 77 or equivalently the 
structural relaxation time r of a supercooled liquid ap- 
proaching its glass transition temperature is called 'frag- 
ile' when it varies according to the so-called Vogel- 
Fulcher-Tammann (VFT) law; and the behavior is called 
'strong' when 77 or t obeys Arrhcnius law |16| . For wa- 
ter, a fragile liquid at room temperature and at mod- 
erately supercooled temperatures, Ito and coworkers 
proposed that a 'fragile-to-strong' transition would occur 
at around 228 K, based on a thermodynamic argument. 
The F-S transition in a molecular liquid like water may 
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FIG. 3: The pressure dependence of the F-S dynamic transi- 
tion temperature, Tl, plotted in the P-T plane (solid circles). 
Also shown are the homogeneous nucleation temperature line, 
denoted as Th , crystallization temperatures of amorphous 
solid water pflj . denoted as Tx, and the temperature of max- 
imum density line, denoted as TMD |2l|. 



FIG. 2: Temperature dependence of (tt) plotted in log((rr)) 
vs T /T or 1/T. Data from ambient pressure, 400, 800, 1200, 
1600, and 2000 bars are shown in panels A, B, C, D, E, and 
F, respectively. 



be interpreted as a variant of kinetic glass transition pre- 
dicted by the ideal MCT where the real structural 
arrest transition is avoided by an activated hoping mech- 
anism below the transition. 

In Fig|2 we report the temperature variation of (tt) 
for water molecules as a function of pressure. It is seen 
that Figs. |2KjI2I3,I2P and[2f) show clearly a transition 
from a VFT law: (r T > = nexp [DT /(T - T )], where D 
is a constant providing the measure of fragility and To, 
the ideal glass transition temperature, to an Arrhenius 
law: (tt) = r\exp[EA/RT], where Ea is the activation 
energy for the relaxation process and R, the gas constant. 
This transition of VFT to Arrhenius behavior, previously 
observed at ambient pressure i s the signature of a 
F-S dynamic transition predicted by Ito et al [17| . In this 
paper, we show the extension of this transition into finite 
pressures. The transition temperature, Tl, as the cross- 
ing point of the VFT law and Arrhenius law, is calcu- 
lated by l/T L = 1/T - {Dk B )/E A . However, in Figs[2E 
andEf , the cusp- like transition becomes rounded off and 
there is no clear-cut way of defining the F-S transition 
temperature. Note that in Fig[2P, there is still a hint of 
fragile behavior at high enough temperature. 

Summarizing all the results, we show in a P-T plane, 
in Fig|31 the observed pressure dependence of Tl and its 
estimated continuation, denoted by a dash line, in the 
pressure region where no clear-cut F-S transition is ob- 



served. One should note that the Tl line has a negative 
slope, parallel to TMD line, indicating a lower density 
liquid on the lower temperature side. This Tl line also 
approximately tracks the Th line, and terminates in the 
upper end when intersecting the Th line at 1600 bar and 
200 K, at which point the character of the dynamic tran- 
sition changes. We shall discuss the significance of this 
point later on. A special feature of the Tl line at the 
lower end should be noted as well. The line essentially 
becomes vertical after around 200 bar and the transition 
temperature approaches a constant value of ~225 K. 

Since Tl determined experimentally is a dynamic tran- 
sition temperature, it is natural to question whether the 
system is in a liquid state on both sides of the Tl, and if 
so, what would the nature of the high-temperature and 
low-temperature liquids be? Sastry and Angell have re- 
cently shown by a MD simulation that at a temperature 
T w 1060 K (at zero pressure), below the freezing point 
1685 K, the supercooled liquid silicon undergoes a first- 
order liquid-liquid phase transition, from a fragile, dense 
liquid to a strong, low-density liquid with nearly tetrahe- 
dral local coordination [22| . Prompted by this finding, we 
may like to relate, in some way, our observed Tl line to 
the liquid-liquid transition line, predicted by MD simula- 
tions of water |j| and speculating on the possible location 
of the low-temperature critical point. 

According to our separate inelastic neutron scatter- 
ing experiments, which measure the librational density 
of states of water contained in 18 A pore size MCM-41- 
S, water remains in disordered liquid state both above 
and below the F-S transition at ambient pressure. Fur- 
thermore, our analysis of the F-S transition for the case 
of ambient pressure indicates that the activation energy 
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barrier for initiating the local structural relaxation is 
Ea = 4.89 Kcal/mol for the low-temperature strong liq- 
uid. Yet, previous inelastic scattering experiments of 
stretch vibrational band of water [23( indicate that the ef- 
fective activation energy of breaking a hydrogen bond at 
258 K (high-temperature fragile liquid) is 3.2 Kcal/mol. 
Therefore, it is reasonable to conclude that the high- 
temperature liquid corresponds to the high-density liq- 
uid (HDL) where the locally tetrahedrally coordinated 
hydrogen bond network is not fully developed, while the 
low-tempcrature liquid corresponds to the low-density 
liquid (LDL) where the more open, locally ice-like hy- 
drogen bond network is fully developed |24|. 

It is appropriate now to address the possible location 
of the second critical point Q ■ Above the critical temper- 
ature Tq and below the critical pressure pc, we expect 
to find a one-phase liquid with a density p, which is con- 
strained to satisfy an equation of state: p = f(p,T). If an 
experiment is done by varying temperature T at a con- 
stant pressure p < pc, p will change from a high-density 
value (corresponding to HDL) at sufficiently high tem- 
perature to a low-density value (corresponding to LDL) 
at sufficiently low temperature. Since the fragile behav- 
ior is associated with HDL and the strong behavior with 
LDL, we should expect to see a clear F-S transition as we 
lower the temperature at this constant p. Therefore, the 
cusp-like F-S transition we observed should then occur 
when we cross the so-called Widom line in the one-phase 
region 25] . On the other hand, if the experiment is per- 
formed in a pressure range p > pc, corresponding to the 
two-phase region and crossing the Liquid-Liquid (L-L) 
coexistence line, the system will be consisting of mixture 
of different proportions of HDL and LDL as one varies 
T. In this latter case, (tt) vs. 1/T plot will not show 
a clear-cut F-S transition (the transition will be washed 



out) because the system is in a mixed state. The above 
picture would then explain the dynamical behavior we 
showed in FigEl In FigEl a clear F-S transition is ob- 
served up to 1200 bar and beyond 1600 bar the transition 
is rounded off. From this observation, the reasonable 
location of the L-L critical point is estimated to be at 
PC = 1600 ± 400 bar and T c = 200 ± 10 K, shown by a 
star in FigEl 

Additionally, in a recent MD simulation using Jagla 
Model Potential by Xu et al .26] , a small peak was found 
in the specific heat Cp when crossing the Widom line at 
a constant p. Meanwhile, Maruyama et al conducted an 
experiment on adiabatic calorimetry of water confined 
within nano-pores of silica gel |27j . It was found that 
water within 30 A pores was well prevented from crys- 
tallization, and also showed a small Cp peak at 227 K at 
ambient pressure. This experimental result further sup- 
ports that the F-S transition we observed at 225 K at 
ambient pressure is caused by the crossing of the Widom 
line in the one-phase region above the critical point [2y| . 
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